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Abstract. The double zeta function was first studied by Euler in response to a letter 
■ from Goldbach in 1742. One of Euler 's results for this function is a decomposition 

formula, which expresses the product of two values of the Riemann zeta function as a 
finite sum of double zeta values involving binomial coefficients. In this note, we establish 
| a q-analog of Euler's decomposition formula. More specifically, we show that Euler's 

decomposition formula can be extended to what might be referred to as a "double g-zeta 
function" in such a way that Euler's formula is recovered in the limit as q tends to 1. 



> 

^ ■ 1. Introduction 

O 

The Riemann zeta function is denned for 9?(s) > 1 by 



OO 

cw==E=- w 



u 

n=l 



n—1 



o 
in 
o 

Accordingly, 

pi ~ i — i 

tf fl '*) : =E^Efct' *( fl )> 1 » »(* + *)> 2, (2) 

^ ■ is known as the double zeta function. The sums (J2J), and more generally those of the form 
C0i,s 2 ,...,s m ) := II^ 7 ' E^( s ^ >n ' ^ = 1 5 2,...,m, (3) 

fel>fc 2 >->fcm>0 j=l 3 j=l 

have attracted increasing attention in recent years; see eg. [21 IH El El El El EE U21 UH1 UH] • 
The survey articles El E21 EH EH] provide an extensive list of references. In (J3J) the 
sum is over all positive integers ki, . . . , k m satisfying the indicated inequalities. Note that 
with positive integer arguments, s\ > 1 is necessary and sufficient for convergence. 

The problem of evaluating sums of the form (J2J) for integers s > 1, t > seems to have 
been first proposed in a letter from Goldbach to Euler jTHj in 1742. (See also ^51 HTj 
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and pH p. 253].) Among other results for (J2j), Euler proved that if s — 1 and t — 1 are 
positive integers, then the decomposition formula 



C(s)((t) = ^( a+ /_ 1 1 )at + a,s-a) + ^( a+ s ^ i 1 )as + a 1 t- a) (4) 

n— n \ / n—r\ \ / 



a=0 v 7 a=0 



holds. A combinatorial proof of Euler's decomposition formula (p£J) based on the Drinfel'd 
integral representations IH 13 El 



<w= n 



dx 4 i dx a 



X 4 J 1 3/ 



l>a:i>— >a; 8 >0 * 1 

£\ _ /" /" TT dx s f tt dx s +t 

J \ _\ x i J 1 ~ x s\ Xi J 1 - X s+ t 

l>xi>->x 8+t >0 1-1 l_s+i 

and the shuffle multiplication rule satisfied by such integrals is given in eq. (10)]. It 
is of course well-known that (jlj) can also be proved algebraically by summing the partial 
fraction decomposition [20, p. 48] ^3 Lemma 3.1] 

1 _y^/q + t-l\ 1 | U + 8 ~ 1\ 1 

a^fc-x)* ^1 i-1 J x s ~ a c t+a ^ \ s-1 ) c s+a (c - xY- a 

K ' a=0 v 7 a=0 V / \ J 

over appropriately chosen integers x and c. (See eg. [H].) 

A g-analog of @ was independently introduced in [TTJ |2TJ 1211 as 



C[s 1 ,s 2 ,...,s m ]:= II r fc .p ' ( ? ) 



where 



fc- 1 - jfc 



Mi-E^T^T' 0<g<l. 

i=o 



Observe that we now have 

C(si,...,s m ) = lim C[si,...,s m ], 

so that (J7J) represents a generalization of (jSJ). In this note, we establish a g-analog of 
Euler's decomposition formula (j3J). 
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2. Main Result 

Our g-analog of Euler's decomposition formula naturally requires only the m — 1 and 
m = 2 cases of (JIJ); specifically the g-analogs of (JTJ) and (j2j) given by 

n (s— l)n n (s—l)n n (k—l)t 

cw=EV- and CM = £ m4* • (8) 

n>0 n>fe>0 L%L«] 9 

We also define, for convenience, the sum 

<p[ 8 ] := V EE_^ = V ^- CM- (9) 

We can now state our main result. 
Theorem 1. If s — 1 and t — 1 are positive integers, then 

s— 1 s— 1— a 



cwcw = E E f°1"* *) (V)a -?) 6 £[*+^- a -6] 

a=0 6=0 ^ ' ^ ' 

+ E E° (° t-T ') C 7 ') (1 - «>' c|s + 4 - ° - 61 



a=0 6=0 

°y (« + 1 - 3 - 1)! (i-„y 

- g (.-j)!( t -j)| -03T)!^ + '-^ 

Observe that the limiting case q — 1 of Theorem ^ reduces to Euler's decomposition 
formula (jlj). 

3. A Differential Identity 

Our proof of Theorem Q relies on the following identity. 

Lemma 1. Let s and t be positive integers, and let x and y be non-zero real numbers. 
Then for all real q, 

JL = V ( a + t - 1 \( t - 1 \ (1 - q) b (l + (q- l)y) a (l + (q- 



+E E 

a=0 6=0 
min(s,t) 



a + s - 1\ /s - 1\ (1 - g) b (l + (q- l)x) a (l + (q - l)y) 1 



E 



s - 1 y V & y y t - a - h {x + y + (q- l)xy) s+a 

a + t-j-l)\ (1 - g) J ' (1 + (? - l)y) s ~ J (l + (g - 



? (s -j)\(t-j)\ (x + y + (g - l)xyY^ 
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Proof. Apply the partial differential operator 



f / dV' 1 1 ( «9 X? 1 



(s-l)!\ dxj (t-l)!\ dy 
to both sides of the identity 

1 1 

-+-+g-l 



xy x + y + (q — l)xy \x y 
Observe that when q — 1, Lemma ^ reduces to the identity 

s-1 / . , -,\ -, t-i 



□ 



1 ^ fa + 1 - 1\ 1 ^-A /a + s - 1\ 1 

x s y t ~ V 1 ~ 1 / ^ s ~ a + y)' +a s ~ 1 / + y) s+ V _a ' 

from which the partial fraction identity (jUJ) (proved by induction in ^J]) trivially follows. 



4. Proof of Theorem [T] 
First, observe that if s > 1 and t > 1, then from (JSJ), 

cwcw = E £ W • W» 

where the inner sum is over all positive integers u and v such that u + v = n. Next, apply 
Lemma^with x = [u] q , y = [v] q , noting that then 

1 + (q — l)x = q u , 1 + (q - l)y = q v , x + y + (q - l)xy = [u + v] q . 

After interchanging the order of summation, there comes 

a=0 fe=0 ^ / V / 



a=0 6=0 
min(s,f) 
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where 



x „(s—l)u„(t—l)v r Jt—l—b)u„av x „(t+a— l)(u+v) a— b— l)u 

or g f a h) — g g g _g _ Or Q __ 



and 



n=l u+v=n L '1 L J< ? n=l M+i>=n L J 5 L Jg 

00 _(t+a-l)n n ~ 1 Js-a-b-l)u 

EQ_ \ ^ <£ 

IriU+a Z-^ L.ls— a— 6 
«=1 ^J<? «=1 L"J<7 

C[£ + a, s — a — b] 



~ q (s-l)u q (t-l)v q (t-j)u q (s-j)v ™ q (s+t-j-l)(u+v) 

i , + J i , \u + V\o 

n=l u+v=n L J y n=l u+v=n L J 

= p[s + t- i]. 



□ 



5. Final Remarks 



In j2^) Zhao also gives a formula for the product C[ S ]CM- However, Zhao's formula is 
considerably more complicated than ours, as it is derived based on the g-shuffle rule [TTj 
satisfied by the Jackson g-integral analogs of the representations ©• Of course, we also 
have the very simple g-stuffle [H] formula CMCM = ([s, t]+([t, s]+([s+t]+(l-q)([s+t-l]. 
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